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R.R.  Put* 

1.  In  the  study  of  the  ocean  and  its  interaction  with  its  environment, 
there  is  a need  for  effective  quantitative  description  of  wave  aotivity  in  the 
ocean.  If  the  surface  is  not  too  irregular,  the  time  history  of  fluctuations 
of  the  water  level  above  some  fixed  point  on  the  bottom  is  a useful  quantity* 
Similar  data  have  been  obtained  for  several  years  by  mean®  of  submerged  pressure- 
sensitive  instruments . The  recorded  time-history  ourve  of  the  subsurfaoe 
pressure  for  a moderate  length  of  time,  twenty  minutes,  for  example,  is  what  wo 
shall  refer  to  as  a wave  record. 

The  questions  we  have  been  concerned  with  aret  What  information  is 
contained  in  a wave  record,  and  how  may  a simm&ry  of  this  information  be  con- 
veniently extracted?  We  shall  describe  a mathematical  model  found  useful 
for  interpreting  wave  records  in  terms  of  probability  distributions  and  energy 
die tributions . Vfe  shall  first  introduce  auxiliary  curves  closely  related  to 
the  original  wave-record  ourve  and  -then  show  how  simple  measurements  on  these 
ourves  can  be  made  to  yield  a description  of  thu  energy  distribution.  We  shall 
then  desoribe  some  statistical  regularities  which  havo  been  observed  in  wave 
records,  and  describe  others  which  have  been  deduoed  as  properties  of  the  speoial 
theoretical  model. 

The  subject  matter  of  this  talk,  as  far  as  observational  data  is 
ooncerned,  is  limited  to  the  time  history  of  pressure  fluctuations  at  a single 
point,  submerged  well  below  the  ocean  surface,  but  in  relatively  shallow  water. 
Some  of  the  methods  and  results  described  are  equally  applioable  to  other  kinds  of 
data. 


2.  In  the  lower  half  of  Figure  1 we  6ee  a segment  of  a typical  wave 
record,  showing  the  pressure  fluctuation  S'  below  the  surface  at  Pt.  Sur, 

California.  The  total  time  interval  show  r ..wits  about  four  and  one-half 
minutes,  within  which  from  twenty  to  twont;  ./e  waves  pass  the  recording  point* 
It  will  be  seen  that  the  distance  the  recording  pen  has  travelled  from  the 
oenter  position,  as  well  as  the  effective  vortical  velocity  and  acceleration  of 
the  pen,  varies  irregularly  from  instant  to  instant. 

Let  us  suppose  that  in  addition  to  the  original  wave-reoord  curve  C0, 
we  oonstruot  a new  tiurve  Ct  representing  the  fluctuating  steepness,  or  slope, 
of  the  ourve  CQ  at  each  point.  This  ourve,  which  represents  the  velooity  of  the 
recording  pen,  may  be  traoed  for  us  by  a machine  if  we  wish.  Applying  suoh  a 
machine  to  the  ourve  C^,  we  shall  obtain  a curve  C2,  representing  the  velooity  of 
the  pen  whioh  traces  C^,  or  the  acceleration  of  the  pen  which  traoes  0o. 

Evidently  this  prooess  of  obtaining  a new  ourve  oan  be  imagined  to  be  repeated 
indefinitely.  We  oall  the  resulting  ourves  the  derived  wave-reoord  ourves. 


* Presented  at  the  1963  Paoific  Section  meeting  of  the  Amerioan  Society  of 
Limnology  and  Oceanography,  Santa  Barbara,  California  17  June  1953. 


TYPICAL  OCEAN  WAVE  RECORO 


Let  us  new  consider  any  varying  curve  and  imagine  a horizontal  straight 
line  drawn  across  it  at  an  arbitrary  height.  Suppose  the  fraction  of  the  time 
the  ourve  spends  below  this  lino  is  measured.  A6  the  line  moves  upward,  this 
fraotion  will  vary,  taking  on  an  increasing  series  of  values  whioh  can  be  plotted 
on  a graph.  A typical  twenty-minute  wave  record  C0  yields  the  graph  shown  in 
Figure  2,  known  ae  its  distribution  curve.  The  vertical  scale  represents  the 
arbitrary  vortioal  level  of  the  line,  the  horizontal  6oale  the  percent  of  the  time 
the  wave-record  ourve  spends  below  that  level.  The  situation  of  the  plotted 
points  in  an  approximate  straight  line  is  characteristic,  and  is  due  to  the 
choice  of  the  horizontal  scale,  vhioh  is  the  familiar  Gauss ian-distribution  scale, 
or  normal -probability  scale. 

The  significant  tiling  about  the  ocean-wave  records  studied  is  the 
apparent  ability  of  the  Gaussian  distribution  scale  to  rectify,  very  nearly, 
the  plots  of  the  distribution  ourves,  not  only  for  CQ,  whioh  has  been  checked 
experimentally,  but  also,  aocording  to  the  theory,  for  G-^  and  C?.  In  each  oase 
there  will  be  some  ohart  level,  oalled  the  mean  level,  below  which  the  curve 
spends  just  fifty  percent  of  its  time.  The  mean  level  for  the  ourves  and 
C2  turns  out  to  be  zero,  while  for  C0  it  oan  be  arbitrarily  set  equal  to  zero 
by  making  all  vertical  measurements  of  the  ourve  height  from  the  mean  level  on 
the  original  wave  record. 

Since  the  mean  levels  are  all  zero,  the  only  feature  distinguishing 
these  straight-line  distribution  ourves  is  their  slope.  We  take  as  the  measure 
of  this  slope  the  difference  between  the  height  of  the  line  at  the  84  percent 
level  and  its  height  at  the  mean  level.  This  difference  is  the  so-called 
standard  deviation  or  dispersion  of  the  distribution.  For  the  wave-record 
ourves  dQ,  C^,  C2,  we  shall  obtain  in  this  way  slopes  denoted  by  <7"0»<yi>o~2» 
respectively. 

3.  Wo  find  an  interpretation  of  these  three  quantities  if  we  next  con- 
sider the  number  of  times,  denoted  by  N0  and  K1#  that  CQ  and  C1#  respectively, 
pass  through  their  mean  levels  (i.e.  through  the  zero  level).  If  the  total 
length  of  the  wave-reoord  is  T seconds,  then  the  average  number  of  zero-level  cros- 
sings per  seoond  will  be  No/T  and  Ni/t,  respectively j multiplying  by  TT , we  de- 
fine the  quantities  a ir  Nq/t  and  U>  x = IT  N]y'T.  Theory  predicts  that  u>Q  =c r \/a~ Q 

and  s 2/<3~l,  relating  the  slopes  of  the  distribution  plots  to  the  mean 
zero-level  orossing  frequencies. 

Let  us  look  further  at  the  ourves  C0  and  C2,  measuring  the  fraotion 
of  the  time  that  the  two  simultaneously  lie  on  opposite  sides  of  their  mean  levels 
--  one  ourve  above  its  mean  level  while  the  other  is  below  its  mean  level,  or 
vioe  versa.  Call  this  fraction  PQ.  If  we  refer  to  a tri gonometrio  table  of 
cosines,  the  value  of  the  quantity  fQ  defined  as  oos  (lTP0)  is,  aocording  to  the 
theory,  a measure  of  the  degree  to  whioh  a point  on  one  of  the  ourves  CQ  or  C2 
determines  the  looation  of  the  corresponding  point  at  the  same  Instant  on  the  other* 
The  quantityf Q is  oalled  the  ooelfioient  of  correlation  between  CQ  and  03. 

Theory  predicts  that  f Q * "(WD/Wi)  = ~(No/Ni)»  lowing  the  negative  nature  of  thie 
correlation.  This  relation  has  been  found  to  check  experimentally  within  a few 
percent,  for  twenty -minute  wave  reoords. 


CUMULATIVE  DISTRIBUTION  FUNCTION  OF  THE  ORDINATES  OF  THE  WAVES 
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FIGURE  3 


4.  A further  interpretation  of  the  distribution  slope  cr  Q for  the  ourve 
C0  is  possible  if  we  examino  the  heights  of  the  peaks  on  this  ourve.  If  these 
heights,  measured  from  the  mean  level,  are  averaged,  the  resulting  quantity^ 
is  Just  one-half  the  average  trough-to-crest  wave  heightyC/^.  We  haveyfr^  ■ M 
(V/2)t  (-  f0)  o~0>  a theoretical  relation  whioh  oheokn  experimentally  to  within 
a few  peroent.  This  is  shown  in  Figure  3,  where./^,  i.e.  2/^,  is  plotted 
against  (Nq/n^)  CT~0  for  each  of  a number  of  twenty-minute  wave  records.  The 
straight  line  representing  the  theoretical  relation  has  elope  (2 Tf)a. 

The  depondenoe  of  the  moan  peak  heighv^y  solely  oncr-0  and  f 
follows  from  the  faot  that  these  two  quantities  determine  the  entire  dis- 
tribution of  peak  heights.  Figure  4 shows  the  theoretical  distribution  of  peak 
heights  for  various  values  of  (*Qu when  is  taken  as  the  vortical  scale  unit 
on  the  reoord  chart.  The  horizontal  scale  gives  the  peroent  of  the  wave  peaks 
■whioh  ocour  below  arv  Jiart  level  M.  This  scale  has  been  ohosen  to  correspond 
to  the  so-oalled  Rayleigh  distribution,  whioh  the  peak  heights  follow  more  and 
more  closely  as  f5  0 approaches  the  value  minus  one.  It  may  be  noted  that  the 
probability  of  a peak  ocourring  below  the  mean  (zero)  level  is  Just  ^(1  + P ), 
which  ia  olose  to  zero  for  numerical  values  of  f0  near  unity. 

Figure  5 shows  a typioal  observed  distribution  of  v/ave  heights,  measured 
from  trough  to  crest  and  platted  against  a horizontal  Rayleigh  scale  on  whioh 
may  be  read  the  peroent  of  v/ave  heights  below  any  height  h.  The  good  fit  of 
the  straight  line  represents  further  agreement,  in  a twenTjy-minute  wave  reoord, 
with  the  theory. 

The  Rayleigh  distribution  enters  into  the  analysis  of  wave  reoords 
in  another  way.  It  con  be  6hovm  that  this  distribution  is  the  one  followed  by 
the  ordinates  on  an  auxiliary  ourve,  associated  with  the  original  wave-reoord 
ourve,  and  known  as  its  envelope.  The  envelope  concept  is  olosely  related 
to  one  aspeot  of  the  air  eady-mentione d theoretical  model. 

5.  The  mathematical  probability  model  approximating  the  wave  reoords 
has  two  aspects.  The  first  is  concerned  with  certain  curve-height,  or  ordinate 
distributions  similar  to  those  already  considered.  The  seoond  is  oonoerned 
with  a theoretioal  distribution  of  wave  energies  over  a continuous  set  of  elemen- 
tary waves  of  all  possible  frequencies. 

The  combination  of  these  two  kinds  of  distributions  oharaoterizes  what 
has  been  oalled  a random  prooess  or  a semi-determinate  function.  The  capaoity 
of  suoh  a theoretioal  model  to  describe  ocean-wave  reoords  is  illustrated  by 
comparing  a sample  taken  from  suoh  a process  and  a typical  wave  reoord.  In 
Figure  1 is  shown  the  two  curves  together;  the  upper  one  was  generated  by  first 
specifying  three  fixed  numbers,  characterizing  the  spectrum,  and  then  con- 
sulting a table  of  random  numbers  published  by  the  CSaubl*i.  Ci  go  University  Press, 

Eaoh  of  the  distributions  appearing  in  tne  general  mathematical  model 
may  be  specified  by  a ourve.  The  ourve  representing  the  energies  of  elementary 
waves  is  called  the  speotral-energy  distribution  curve.  This  distribution  of 
energy  is  of  fundamental  iraportanoe  for  the  study  of  the  generation,  propagation 
and  meohanioal  effeots  of  waves.  Its  direct  computation  presents  oortain  prao- 
tioal  difficulties  and  would  appear  to  be  too  time-oonsuming  for  routine  wavo- 
date.  analysis.  Our  problem  is  the  partial  characterization  of  the  speotral- 
distribution  ourve  by  means  of  a very  few  selected  quantities. 


RELATIVE  DENSITY  OF  DISTRIBUTION  OF  POWER  IN  SPECTRUM 
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In  Figure  6 are  shown  certain  members  of  a mathematical  family  of 
speotral  distribution  curves  whioh  are  completely  characterized  by  four  numbers 
each*  Three  of  these  numbers  are  the  same  for  each  of  the  speotral  distribution 
curves  shown.  These  common  parameters  are;  (1)  the  total  area  enolosed  by 
the  curve,  whioh  turns  out  to  b e o~  (2)  the  mean  speotral  frequency, * 

averaged  over  the  speotral  energy  d?stribution,  and  (3 ) the  standard  deviation 
or  bandwidth  o~ ^ of  the  spectral  energy'  distribution.  The  fourth  number, 
denoted  "by <£3,  is,  for  the  family  in  Figure  6,  a measure  of  the  degree  of  a- 
symmetry  or  skewness  of  the  speotral-distribution  curve.  It  is  seen  that  there 
is  a marked  degree  of  asymmetry  whenfifc-g  = 2.0. 

It  is  clear  that  the  higher  the  speotral  energy  distribution  curve, 
the  greater  the  enolosed  area,  so  that  <r*  Q 2 represents  the  total  energy  for 
the  wave  record.  Likewise,  the  spectral  bandwidth  is  a measure  of  how 
widely  the  energy  is  dispersed  over  all  possible  frequencies.  A.  quantity  of 
some  interest  is  the  relative  bandwidth  £ , defined  as  the  ratio 
measuring  the  degree  of  spread  of  the  energy  relative  to  the  midband  location 
at  which  it  is  centered. 

6.  It  is  possible  that  for  oertain  purposes  a spectral  distribution  may 
be  adequately  described  by  means  of  a few  quantities,  provided  something  is 
known  about  the  general  shape  of  the  spectral  curve.  With  present  knowledge 
of  typical  spectrum  shapes,  the  quantities  just  described  would  appear  to  be 
useful  for  a general  oharaoterii-ubion  of  the  energy  distribution.  In  fact, 
for  simple  speotra  the  three  numbers  <T“02,  , and  , may  prove  to  be 

adequate  for  many  purposes  if  supplemented  by  one  or  two  other  parameters  (such 
as  cC  3)  to  indicate  speotral  shape. 

We  can  now  describe  how  this  kind  of  information  about  the  speotral 
distribution  may  be  obtained  from  the  wave  record.  This  doseription  will  re- 
quire the  concept  of  the  envelope  curve  mentioned  earlier?  Y/e  ehall  show  later 
how  in  special  cases  approximate  information  about  the  spectral  distribution 
may  be  obtained  by  considering  only  the  previously-defined  derived  wave-reoord 
curves. 

For  a regular  wave  system  vshose  speotral  energy  is  concentrated  at 
a single  frequency,  one’s  intuition  rightly  suggests  as  the  envelope  curve, 
the  horizontal  straight  line  passing  through  all  the  peaks  on  the  wave-record, 
which  has  in  this  case  a simple  sinusoidal  form.  For  more  general  wave  systems, 
the  envelope  ourve  may  be  determined  mathematically,  but  its  exact  location 
on  the  chart  is  not  as  easy  to  visualize.  It  may  be  said,  however,  that  if  the 
wave-reoord  curve  is  not  too  irregular,  the  envelope  curve  oomes  close  to  each 
peak. 

Suppose  for  a given  wave-reoord  curve  the  envelope  is  drawn  and  a 
new  curve  then  constructed  at  each  point  in  the  following  way;  Form  the 
squares  of  the  heights  of  the  ourve  and  of  its  envelope  at  each  point,  subtraot 
the  former  from  the  latter,  and  extract  the  square  root  of  the  difference.  The 
result  is  the  numerical  value  of  the  height  of  a point  on  a now  curve  said  to 
be  conjugate  to  the  original  ourve.  The  entire  operation  leading  to  the 
drawing  of  the  conjugate  curve  may,  in  certain  cases,  be  carried  out  by  ma- 
chines with  sufficient  aocuraoy. 


7.  If  C0 * and  Ci*  a r«  the  ourvea  conjugate  to  CG  and  Ci  respectively, 
these  flour  ourves  may  be  paired  together  in  six  different  ways.  If  for  eaoh 
pair  we  measure  the  fraction  of  the  time  the  two  curves  lie  on  opposite  sides 
of  their  mean  levels,  we  find  that  except  for  the  pairs  (Cq,C^*)  a£U*  (C.,  C Q*)» 
this  fraotion  is  always  one-half,  independent  of  the  speotrum  of  the  original 
curve  C0.  For  the  remaining  two  pairs,  this  fraotion  depends  on  the  speotrum 
of  C0#  containing  in  both  oases  essentially  the  same  speotral  information, 
whioh,  it  turns  out,  is  just  what  is  needed  for  c lputing  the  moan  spsotral 
frequency,/4^  for  the  • ive-record  ourve  CQ. 

l^t  PQ*  be  the  fraotion  of  the  time  the  ov'-ves  0o  and  C^*  lie  on 
opposite  sides  of  their  mean  levels*  Then  *,  the  value  given  in  a table 
of  oosinee  for  the  angle  (Tt^P  *),  will  be  the  coefficient  of  correlation  be- 
tween CQ  and  Cj*.  This  quantity,  together  with  the  mean  zero-level  crossing 
frequency tt>Q , determines  the  mean  speotral  frequency*  In  faot,  wo  have 
^ p0*n-  Further,  it  may  be  shown  thatx^3“(°0  whore  to  ^as 

before,  is  the  mean  wave-crest  frequency,  measured  by  counting  all  peaks  ao 
matter  how  high.  The  correction  factor  (*f  o fo  *)  nooessary  to  oonvert  from 
mean  poak-to-peak  frequency,  CO^t  to  mean  speotral- distribution  frequenoy,/^., 
nay  be  as  small  as  0.7  for  a typical  pressure  record.  The  quantity  ffo  fo  ) 
is  in  all  practical  cases  less  than  unity  and  smaller  that  the  similar  quan- 
titv  f Q*,  whioh  is  the  correction  factor  required  when  uj0  is  used. 

It  will  be  seen  that  if^/^^is  regarded  as  a suitable  measure  of 
mean  wave  frequency,  then toQ  and  co-^  probably  would  not  be,  since  eaoh  of 
these  two  quantities  includes  information  about  the  speotral  bandwidth  aa  well 
as  about  the  center  position  of  the  speotral  band.  The  extent  to  which  fre- 
quency bandwidth  thus  enters  into  the  mean  ^ero-crossing  frequency  («o). 

whioh  might  otherwise  appear  to  be  a simple  measure  of  midband  frequency,  depends 
upon  the  relative  bandwidth  (Siti)  the  spectral  distribution. 

» 

8.  If  we  consult  the  trigonometric  tables  once  more,  finding  the  numer- 

ical value  of  the  tangent  of  the  angle  ( 7f  PQ*),  we  obtain  the  relative  band- 
width itself.  The  smaller  this  tangent,  the  mere  nearly  the  wave-record  curve 
will  lc  a single-frequency  wave  system.  The  smaller  this  quantity,  the 

more  closely  correlated  will  be  all  the  derived  and  conjugated  curves  ob- 
tained from  CQ  — with  the  exception  of  those  which  simultaneously  lie  on  op- 
posite sides  of  their  means  exactly  one  half  of  the  time,  regardless  of  C0'a 
speotral  distribution.  Figure  7 shows  how,  for  a rather  wide  class  of  spectral 
distributions  including  thote  shown  in  Figure  6,  the  values  of  f Q and  Sat 
determine  one  another  with  inoreasing  preoisicn  as  the  relative  bandwidth  de- 
creases. The  region  between  the  two  curves  shovm  represents  the  possible 
variation  of  these  quantities,  the  upper  ourve  corresponding  to  e£.  3 * 0.0,  the 
lower  to  ct  3 a +2.0.  As  the  relative  bandwidth  doorcases,  the  distribution 

of  peak  heights  becomes  mere  nearly  like  that  of  the  envelope.  The  envelope 
ourve  itself  then  approaches  the  wave-reoord  curve  at  each  peak  of  the  latter. 

The  oase  of  small  relative  bandwidth  is  significant  since  most  modera- 
tely deep  pressure  records  seem  to  fall  into  this  category  with  sufficient 
precision  to  allov/  certain  narrow-band  approximations  to  be  made,  in  most 
cases  with  an  error  of  a few  percent.  V/e  then  have  the  following  approximate 
formulas  requiring  no  knowledge  of  the  conjugate  wave-record  ourve.  In  terms 
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9.  The  random  process  interpretation  of  wave  reoords  deals  directly  with 
the  properties  of  the  unobserved  part  of  a wave  system,  for  ifc  makes  the  baBlo 
assumption  that  the  wave  reoord  is  only  a sample  from  a oontinuing  prooess. 

This  approach  has  been  used  to  prediot,  by  a least-squares  method,  the  behavior 
of  waves  at  one  time  and  plaoe  from  the  behavior  of  waves  observed  at  a time 
and  place  nearby.  In  the  lower  half  of  Figure  8 are  shown  the  wave  reoord  ob- 
served, and  that  predicted, for  the  same  point,  using  the  wave  reoord  observed  at 
a point  approximately  200  feet  seaward.  This  prediction  was  made  after  an 
Initial  analysis,  not  of  these  two  \wave  reoords,  but  rather  of  two  wave  reoords 
made  at  these  points  thirty  minutes  earlier.  The  prediction  method  requires  the 
determination  of  but  four  numbers  from  eaoh  of  the  speotral  distributions 
involved.  The  upper  half  of  Figure  8 shows  the  results  when  the  four  numbers 
are  determined  from  these  wave  reoords  themselves,  rather  than  from  earlier 
wave  reoords  at  the  same  points. 

10.  A word  of  qualification  should  be  added  to  the  theoretical  results 
that  have  been  indicated.  All  of  the  quantities  whose  relations  to  eaoh 
other  have  been  described  are  assumed  to  have  been  oaloulatsd  from  wave 
reoords  of  suffioient  length  to  minimise  random  sampling  errors.  There  is  some 
evidence  that,  for  describing  ooean  waves,  a twenty-minute  reoord  may  be  ade- 
quate, although  for  some  purposes,  longer  reoords  may  be  neoeBsary.  Shorter 
wave  reoords,  of  course,  reduce  the  likelihood  of  significant  changes  in  the 
wave  spectrum  during  the  ti.aa  oovered  by  the  reoord,  due  for  example,  to 
meteorologioal  conditions. 

It  may  be  suggested  that  the  quantities  described,  when  continuously 
obtained  from  field  installations  by  routine  mechanised  analyBiB,  will  be  a 
useful  supplement  to  the  results  of  an  oooasional  more  refined,  time-consuming 
analysis  of  the  wave  spootrum.  Quantitative  information  would  then  be  availa- 
bel  during  the  life  history  of  a storm  describing  changes  in  the  wave  speotr'"* 
as  they  ooour.  Numerical  results  would  be  available  for  investigating  both 
the  generation  of  waves  and  their  ultimate  effect  on  the  coast  and  other  struc- 
tures, as  well  as  for  the  establishment  of  wave-olimate  conditions  at  various 
looations  and  Beasons  of  the  year. 
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